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Abstract 

In this work, a chain of A/ identical two-level atoms coupled with a quantized electromagnetic field, initially prepared 
via a single-photon Fock state, is investigated. The A/-particle state amplitude of the system is calculated for several 
space configurations of the atoms in the Weisskopf-Wigner approximation. It was shown that the space configuration 
of an atomic chain, the total number of atoms, and even the available volume for the field modes define the behavior 
of the system state amplitude with time. Applying the condition of 'cyclic bonds', presented in this work, to the 
elaborated theory allows to describe the system time evolution, practically, for any space configuration. 
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Background 

The collective absorption (emission) of photons by an 
ensemble of identical atoms provides valuable insights 
into the many-body physics of photons and atoms' 
(quoted from [1]). Taking into account the quantization 
of electromagnetic field, many fundamental and interest- 
ing properties of the coupled systems of atoms and field 
are revealed. For example, when the average distances 
between atoms are much less than the resonant transi- 
tion wavelength of emitted (absorbed) light, the cooper- 
ative coupling leads to a substantial radiative shift of the 
transition energy and significant change in decay rate of 
the ensemble state. The latter was revealed through the 
various theoretical (for example, some relatively modern 
researches in [2-5]) and experimental investigations (see 
starting, for example, from [6,7] to the modern applica- 
tions like described in [8] and impressively effective exper- 
imental realizations as in [1]). Some peculiar behavior in 
spontaneous emission is proper even in a system of atoms 
which can have a relative distance larger than the emission 
wavelength (see, for instance, [9]), and initially, only one 
atom or one-photon state is excited, as discovered in this 
paper. 
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In the present paper, a system (chain) of N identical 
two-level non-interacting atoms, prepared Via a single- 
photon Fock state' in the one- or two-mode resonator, 
is investigated. The main goal of the paper is to obtain 
the information about the state of electromagnetic field 
and atomic system (chain) in a Weisskopf-Wigner approx- 
imation (see [10] chapter 6, page 206 and some com- 
ments in [11]). The calculations of the state amplitudes 
of the atomic system are made for several approxima- 
tions in resonator (cavity) characteristics and for sev- 
eral types of space configurations of the atoms in the 
chain. 

In this work, we study the case, in which the distances 
between atoms are quite large, so that the average dis- 
tances between atoms are greater or in the same order 
than the 'resonant transition wavelength. Therefore, we 
prepare an ensemble of N two-level atoms initially in 
ground state, and a single mode of the radiation field is 
excited in a 'Fock' state (so called one-photon state). This 
is the case of a purely monochromatic wave with zero 
line width under the consideration. A laser output in sin- 
gle mode operation can approximate this situation due 
to its high degree of monochromaticity (small line width) 
for instance. The mode of electromagnetic field is speci- 
fied completely by giving its wave vectors ko with atomic 
transition frequency co = c|kol and its polarization ; 
(/• = 1, 2). 
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The main feature, differentiating our research from 
others in this domain, is the developed direct and con- 
sistent solution to the A/^-particle equations, describing 
the time evolution of the N atomic probability state 
amplitudes. Besides, in certain sense, we explained the 
nature of the widely used Weisskopf-Wigner approxima- 
tion that was not found in the reviewed by us scientific 
literature. 

The goal of this paper can be formulated as an attempt 
to propose an adapted and simple in practical use the- 
ory, for example in the highly applied nanoscale physics. 
The proposed theoretical material requires correspond- 
ing experimental verification. As an idea of an application, 
the model system can be realized on atomic (develop- 
ing the method proposed in [1] for the nuclei of ^^Fe 
in certain composites, but this time for a visible region), 
chains of trapped ions (like in [8]), and molecular struc- 
tures for further developing such techniques like FRET 
(described for instance in [12]), atomic chains like carbyne 
loops (for example, [13]), and microhole array synthe- 
sized by femtosecond laser radiation (see [14], for an 
instance). 

Let us first provide below some general theoretical 
premises. More detailed derivations of the corresponding 
mathematical model can be found in [11]. 



and EM field is in Fock state 1 1^^ ) (that presents one pho- 
ton with the wave vector ko), we look for a solution of 
the corresponding Schrodinger equation in the interacting 
picture in the following form: 

N 

=Y^fia{t)\bib2...aa...bN0) 

a=l (1) 

+ XI yk.j(i)\hb2 • • • ^Tvlk,;) 

with the initial conditions: 

Pa(0) = 0, ]^,;(0) = 5k,ko> (2) 

where Sy^^y^Q is Kroneker s delta symbol. Sy^^i^^ = 1 if k = 
ko, and (5k,kQ = 0 if k 7^ ko- fia(t) (oi = h^N) and yi,j(t) 
(j = 1, 2) are the ath. atom excited state amplitude with 
the others in the ground states and excited Fock field state 
amplitude of the yth polarization with the wave vector k, 
accordingly. 

Then, the corresponding Schrodinger equation in 
the interacting picture yields the following system of 
equations: 

^ot{t) = /X^*(k,;)]/k,y(0exp(-/(v^-6t))t+/kro,); (3) 



Methods 

The equations of motion for the state amplitudes 

We have assumed that the atomic energy levels have 
no linewidth, so that, only if a; = co/^q = A£^^//i, 
the atoms can be able to absorb a photon. Obviously, 
this is an unrealistic case since it is impossible to have 
a completely monochromatic wave. In addition, for the 
case of the Fock initial state, in which we measured the 
energy precisely of the mode, the average electric field 
will be zero. In the forth of the law of energy conser- 
vation, an emitted photon will correspond to the same 
frequency co/^^ (we can say it will occur with a high prob- 
ability after a quite long time interval if the system has a 
damping). 

Therefore, consider a collection of N identical atoms, 
at positions ri, . . . , Tq,, . . . , r^, coupled to a one mode 
electromagnetic (EM) field. Each atom a = 1..N is 
assumed to have only the two states \a)a and \b)a, sep- 
arated by energy = E^a — E^a = hco. In the dipole 
approximation, the closed conservative system of identi- 
cal atoms with the electromagnetic field in a cavity can 
be described by the Hamiltonian consisting of free atoms 
and electromagnetic field items with dipole-field cou- 
pling between the atoms and the electromagnetic field 
modes. 

Inasmuch as at the initial time moment t = 0 all atoms 
a = 1,.N of the ensemble are in the ground state \b)a 



N 

= ^ ^g8(^J)P8(t) exp(/(v^ - a))t - /kr^), (4) 

where 



where 

Pa -ek,; = I Pa I COS ^k,;. (6) 

Here, ^k,/ is the angle between dipole transition vector 
Pa (more accurately, non-diagonal dipole matrix element) 
and the yth unit polarization vector ek,y (j = 1,2 and 
ek,y • k = 0). y is the available by the system of atoms and 
field space volume. The frequencies v/^ correspond to the 
modes with the module of the wave vectors k equal to |k|. 

Therefore, substituting Equation 4 into 3 and differenti- 
ating one more time, after applying the Weisskopf-Wigner 
approximation (details in [11]), we can derive the follow- 
ing system of evolution equations: 

^2 ^ ^ 

^Pa(t) = -J2Ps(t)^cc8 -2Da-Pa(tl (7) 

where 

^0.5= ^«(M^5(k,;)exp[/k(rc -r^)]. (8) 

;,|kN/:o 
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And the decay rates in the approximation can be 
estimated by the formula: 



1 Stt 27t 



2 3 2hsQ \27zc 



1 



IpI^^^. Bc = Y^ Pait) cos(kr,,) = C (exp (^^2+^) - exp {Q2-t)) , 

a 

(15) 



2 Stt Usqc^ 



(9) 



The coefficient Da {oi = 1..N) describes the respective 
rate of decay for ath atom excited state. Note, that the 
'non-resonant' items for the particle with distinguished 
from a indexes were disregarded in here in an assumption 
of quite large interatomic distances (see details in [11]). 

Results and discussion 

An atomic chain with cyclically distanced atoms: 

N 

J2 sin(kroj) cos(kroj) = 0 

a 

Next, we try to make the calculations, using here the par- 
ticular case of space configuration for the system atoms 
field. Below, for simplicity, only one polarized mode (j = 
1) of the resonant field modes is taken into account with 
the common parameters ^q; and Pa for a = 1,.N: 



gaO^J) = ga=g>0 



and 



Pa -ek,; = IpI 



(10) 



(11) 



for |k| = ko.ln other words, the space angle distribution 
for the components is disregarded here, assuming the 
direction of the transition dipole moment for any atom 
in the system coincides with the photon polarization in 
absorbing or emitting a resonant photon. 

Then, from the system of Equation 7, in the case of a 
cavity with two resonant modes k = ibko and identical 
atoms with = D for a = l..A/^, one derives that 



N 



Pait) = -2/ J2 ^5^^) cos(k(r« - rs))-2D-Ja(t)^ 

8=1 



Using the notation 



N 



Beit) = ^Pa{t) COS {VXa) 



(13) 



N 



and the 'cyclic' condition ^ sin(kro;) cosCkr^,) = 0 yields 
the following relatively simple linear differential equation: 



^2 ^ 

Beit) = -2/ ^cos2(kr,,)5,(0 - 2D-Bc(t). 



dfi 



(14) 



Therefore, taking into account the initial conditions 
yg^(O) = 0 for a = 1..N, the solution of the above equation 
is as follows: 



where ^22=^ X! cos'^(kro;) and 



(16) 



Q.2± = -D± ^D^ - Qj. 
By analogy, 

Bs = J2 sin(kr,,) = C (exp (q2^^ - exp (^2-^)) ^ 

a 

(17) 



where = ^ 2J2 sin'^(krQ;) and 



2± = -D± yjn'^ - 



(18) 



It is easy to see, that ^ + ^2 — 2^^A/^. 

The field probability amplitudes can be obtained using 
the subsystem of Equation 4 of the full 'conservative' sys- 
tem of Equations 3 and 4. Therefore, substituting (15) and 
(17) into the Equation 4, and then taking into account the 
restrictions ^qj(O) = 0 for a = 1..N, we obtain that 

KkCO = 2ig [C {^^2+/ (^^2+, t) - Q2-f (^2-, t)} 
-id [q^-/ (^2+^ - ^2 J (^2-^ ]] + 1> 



(19) 



and 



y_k(^) = 2ig[c {Q2^f (Q2^a) - Q2-f(^2->t)} 



(12) where 



exp(^) ^ /Qt 



(20) 



(21) 



Note, here, we neglected the possible space angle dis- 
tribution for the direction of the resonant wave vector 
k. 

Inasmuch as cos (k(rc^ — ^s)) = cos (kr^;) cos (kr^) + 
sin (kro;) sin (kr^), then, after substitution of the found 
superpositions (15) and (17) into the initial Equation 12, 
we derive the following integrable differential equation: 



d 



^Mt) + ^D-Mt) = - 2g' {cos (kia)Bc(t) ^22) 
+ sin(kroi)Bs{t)}. 
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Integrating the left and right sides of the equation above 
(22) over time yields 



at 
where 

Ta it) = - 2i {cos (kr„) FciO + sin (kr„) 
+ ^/6«(0) + 2D;6„(0), 



(23) 



(24) 



and 



(25) 



According to the definition of the functions Fc,s{t) 



Fcit) = ^ { ^ [e^'P - 1] - ^ [e^'P - 1] j ; 

(26) 



and 

Fs(t) = C 



(27) 



The solution of such linear first order differential 
equation, like (23), has the form: 

fia (t) = ^^^^^2Di) / ^""^ "^^^ ^^^^ 

The integration in the last expression can be performed, 
yielding 

= - 2/ Icos (krj C I ^ r ^ (^((2i)+^2+)^) _ i) 

_ J_ (^^(2^0 _ 1^1 _ J_ r ^ /g((2D+^22-)0 _ ^ 

2D V /J ^^2- L^2- +2D V / 

-^(.(^^^)-l)]j+sin(krJC 



^^2+ + 2D 



^ / ((2D+4_)^) _ \ 

+2D\ J 



+ (^^a(0) + 2^(0)) ^ - i) + Co. 



Therefore, 

= -2^2 j COS C [// (^22+, A 0 - ^ {^2-,D, m 

+ sin (kra) C' |^// (^i^'j.^, D,tj-H (q'2_,D, tjj | 

+ (|;S„(0) + 2Dfia{0)) ^ (1 - + Coe-(2^*>, 

(30) 

where 



H(Q,D,t) = 

^2 



( 



^2 + 2Z) 
2D V / 



(31) 



The initial condition fiai^) = 0, for a = I..N, sets 
the coefficient Co equals 0. The initial time derivative 
^ Pa (0) can be determined, for example, if the system of 
Equation 3 from the initial 'conservative' full system of 
Equations 3 and 4 is chosen as a basis at the time moment 
^ = 0. Then, the initial condition for the field state ampli- 
tude Yis^iO) = 1, where k = ko, sets the time derivative 
^ Pa (0) to the following expression: 

d 

(32) 



= ig {cos (kr^;) + / sin (kr^j)} . 



Now, the question arises how to choose correctly the coef- 
ficients C and C . First of all, the choice has to satisfy 
the limitations on the probability amplitude, yielding the 
corresponding probability limited above by unit (the sum 
of all the modules squared of the introduced amplitudes 
equals unit probability). Secondly, the solution with the 
coefficients have to be consistent with the model decay 
(damping). 

We observe that, formally, when the real part of the 
variable ^2 is a negative quantity, that is Re(Q) < 0, the 
introduced functions H and/ have the following limits for 
quite long time intervals: 



lim H(Q,D,t) = — 

t^oo 2DQ 



, when Re (^2) < 0; (33) 



lim / (Q, t) = ^, when Re (Q) < 0. (34) 

t^oo 2Q^ 



Then, 



lim Kk(0 = ig\c\-^- 

t^oo I [Q2- ^2+ 



iC 







{"2- 


'^2+[ 



■I' 



1 


1 




^^2- 


^2+ ^ 





(29) 



+ 1; 
(35) 



(36) 
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^0 



sin (kva) 



(37) 



As for an open system, in our case, it should be expected 
for a quite long time interval the total electromagnetic 
energy of the atoms-field system to be emitted into the 
subsystem causing the state damping. Therefore, let us 
define the coefficients C and C in the following manner: 



C = 



and 



C = - 



^2-^2+ 



2^(n2+-«2-)' 



(38) 



(39) 



Then, after substitution into the expressions for the time 
limits, one derive the logical finale of the system evolution: 



lim Yi^it) = 0; 



lim K_k(0 = 0; 



lim ^0.(0 =0. 



(40) 



(41) 



(42) 



The possible space configurations of the atomic system, 
satisfying the condition of circularity', can be easily found. 
For example, the set s3al (the notation s3al' is just intro- 
duced here): kri = k • ri = kr2 = and krs = tt. 
As an instance, it can also be the set s3a2: kri = 
kr2 = and krs = + 27r. A space configuration for 
five atoms can be represented, for example, by such set like 
s5al: kri = k • ri = ^, kr2 = ^, krs = ^, kr4 = ^, 
and krs = Specifically, a combination, such like 

kvi = ± ^ -\- Oi (^) + 2i7T with / = h,N, n > 2, and 
|o/ (^)l << ^, can generate the necessary magnitudes 
of the characteristic system frequencies ^^2 and (that, 
actually, are the corresponding Rabi frequencies), compa- 
rable with the given magnitude of the decay coefficient 
D. 

Below we depict the atomic system behavior in the sev- 
eral introduced above configurations. Note, that the cited 
thereby Rabi frequencies were calculated in the SI sys- 
tem of units with the following notations: h ^ 1.05457 x 
10"^"^ J sec/rad; the electric permittivity of free space 6o ^ 
8.8542 X 10~^^ F/m; the speed of light in free space c = 
299792458 m/sec; resonant wavelength close to the D2- 
line of a sodium atom Xd ^ 589.29 x 10~^ m; correspond- 
ing circular (in radians per second) resonant frequency 
a),es = ^ 0.101747 x IO^^tt rad/sec; non-diagonal 
so called 'transition dipole matrix element (in the same 
order as for the D2-line transition, that is about 1 Debye) 
Pex = 1 X 3.33564 X 10"^^ C m. For instance, if the avail- 
able for the system of atoms and field volume has the 



value equal to V = 0.001 m^, then ^ = Pex^ ^ 

77.8597 V^r rad/sec. 

Assume, for example, the available volume V = 
10~^^ m^ is somehow filled by the set s3al with D ^ 
10^ rad/sec, initially coupled with one-photon Fock state. 

Then, g = Pex^j^^ ^ 77.8597 x 10^ 0r rad/sec, 
^^2 ^ 155.7195 X 10^ rad/sec, and Q2 ^ 77.8597 x 



10^ V^TT rad/sec. The corresponding graphs for probabil- 
ity to find each atom in the excited state are shown in 
Figure 1. 

Let us see what happens when the available volume is 
increased by one order. This yields V = 10~^^ m^ with the 
same three atoms {D ^ 10^ rad/sec) of the configuration 
s3al. Then,^ ^ 24.6214x 10^ 0r rad/sec; ^^2 ^ 49.2428x 
10^ rad/sec and ^ 24.6214 x 10^ V^tt rad/sec . 
The corresponding graphs for each atom excited state 
probability are depicted in Figure 2. 

Suppose now that the available volume isV = 10~^^ m^, 
somehow filled by the set s5al with D ^ 10^ rad/sec 
initially coupled with one-photon Fock state. Then, g = 

Pex^f^^ ^ 77.8597 X 10^ Vtt rad/sec; ^^2 ^ 77.8597 x 

10^ rad/sec, and ^^2 ^ 155.7195 x 10^ V^tt rad/sec. 
The corresponding graphs for each atom excited state 
probability are shown in Figure 3. 

And again, let us see what happens when the avail- 
able volume is increased by one order. This yields V = 
10~^^m^ with the same five atoms {D ^ 10^ rad/sec) 
of the configuration s5al. Then, g ^ 24.6214 x 
10^ rad/sec; ^^2 ^ 24.6214 x 10^ V^tt rad/sec, and 



^^2 ^ 49.2428 X 10^ V27r rad/sec. The corresponding 
graph is in Figure 4. Note that the graphs Figures 3 and 
4 of excited state probabilities are for the chosen three 



Probability of an atomic excited 
state as a function of time (in sec. 




2e-07 2.5e-07 3e-07 



space phase = Pi/6 

space phase = 2Pi/3 
space phase = Pi 

Figure 1 Time evolution of \Pci(t)\^- V = 10"^^ m^. Atoms are 
arranged in the set s3al with D ^ 1 0^ rad/sec. The bold solid line 
represents the atom with the space phase kri = 7^/6, the dot line is 
for the space phase Wxj = In/S, and the thin solid line corresponds 
to kr3 = TT. 
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Probability of an atomic excited 
state as a function of time (in sec.) 




t 



space phase = 2Pi/3 

space phase = 19Pi/6 

space phase = 5Pi/2 

Figure 4 Probability \Pait)\^- ^ = 10~^^ m^. Atoms are arranged in 
the set s5al with D ^ 1 0^ rad/sec. The bold solid line represents the 
atom with the space phase kri = 27r/3, the dot line is for the space 
phase krs = 1 97r/6, and the thin solid line corresponds to 
krs = 571/2. 





Probability of an atomic excited 
state as a function of time (in sec.) 


0.03 r 




0.025 T 




0.02 - 




0.015 r 




0.01 T 




0.005 T 




0 ^ 


2e-07 4e-07 6e-07 8e-07 le-06 

t 

space phase - Pi/6 

space phase - 2Pi/3 

space phase - Pi 


Figure 2 Atom excited state probability \Pcc(t)\'^- V = 10"^^ m^. 


Atoms are arranged in the set s3al with D ^ 1 0^ rad/sec. The bold 


solid line represents the atom with the space phase kri = 7r/6, the 


dot line is for the space phase kr2 = 27r/3, and the thin solid line 


corresponds to krs = tt. 



atoms with the following phases: kri = krs = and 
krs = '-f. 

As it was supposed in the derivative of the differential 
equations with the damping items such like (12) (see the 
details in the work [11], the available volume V for the 
system of atoms and field defines the available' modes 
for the electromagnetic field. The value of volume V can 
determine one of the inequalities D < Q2 and D > Q2 
(D < Q2 ^ > ^2)' therefore defining the char- 
acter of the system relaxation. Such fundamental system 
property was illustrated in the figures. It is interesting 
to note that increasing the system volume V, therefore 
increases the available' number of quantized field modes, 
the maximum probability to find an atom in its excited 
state decreases. Other interesting feature, shown in the 
proposed graphs, is the different character of relaxation 



Probability of an atomic excited 
state as a function of time (in sec.) 




t 



space phase = 2Pi/3 

space phase = 19Pi/6 

space phase = 5Pi/2 

Figure 3 Atomic excitation probability \fiait)\^ as a function of 
time. V = 1 0~^^ m^. Atoms are arranged in the set s5al with 
D ^ 1 0^ rad/sec. The bold solid line represents the atom with the 
space phase kri = 27r/3, the dot line is for the space phase 
krs = 1 971/6, and the thin solid line corresponds to krs = Stt/I. 



for each excited atom. The latter depends, as shown here, 
on the space phase kr^;, where a = 1„N, 

On this note, therefore, let our narration to come to the 
following conclusions, in short. 

Conclusions 

Thus, in this work, we investigated a chain of N identical 
two-level long distanced atoms prepared Via a single- 
photon Fock state! The functional dependence of the 
atomic state amplitudes on a space configuration and time 
is derived in the Weiskopf-Wiegner approximation. 

It was shown that in increasing the system volume V, 
the maximum value of probability to find an atom in its 
excited state decreases. The feature can be experimentally 
investigated at the proposed nanoscale limit for the space 
configuration of atoms. 

Hence, the Weiskopf-Wiegner approximation was 
revealed through the provided application to the many- 
body system at the nanoscale limit for the atomic space 
phases. The found solution (30) cannot be counted as a 
particular one, or as a limit of such, for the initial systems 
of Equations 3 and 4 that represent only a closed con- 
servative system of atoms and an electromagnetic field. 
Thus, we can say that the model described in this work , 
besides the atoms and the electromagnetic field, implicitly 
contains a third participant guaranteeing a total system 
relaxation with time. It is interesting to note here that the 
complete' decay of the system excitations was strongly 
imposed by the choice of the coefficients C (38) and 
d (39). 

The methods, described in this work, of solving the 
system of linear differential equations can be applied 
even for more general situations when the boundary 
circular' conditions are not satisfied. In certain cases, 
the problem allows to extend the system by adding a 
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subsystem of non-sufficient number of atoms (or an 
atom when the expression sin(krc,,) cos(krc,,) | < 1 

is actual) without influencing the main characteristics 
under our interest. Besides, it is interesting to investigate 
the dependence of value of the mentioned construction 
sin(krQ,) cos(krQ,)| on the location of the origin of a 
coordinate system. 
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